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Abstract 

^ ' We show that the connected-moments polynomial approach pro- 

^_S ■ posed recently is equivalent to the well known Rayleigh-Ritz variation 

■tH- . method in the Krylov space. We compare the latter with one of the 

original connected-moments methods by means of a numerical test on 
t~^ ■ an anharmonic oscillator. 

o 

00 

• • '. 1 Introduction 

> ■ 

/\ • The t-expansion [1] has motivated a great interest in the application of the 

j^ ■ connected-moments expansion (CMX) [2] and its variants [3-5] to quantum- 

mechanical models. In spite of its limitations [3-8] the CMX has proved 
useful for the study of many physical problems [2,3,6,7,9-12]. 

Recently, Bartashevich [13] proposed the connected-moments polynomial 
approach (CMPA) that yields approximate eigenvalues to all states as roots 
of a simple polynomial function of the energy with coefficients that depend 
on the moments of the Hamiltonian operator. 

The purpose of this paper is to show that the CMPA is simply the 
Rayleigh-Ritz variation method (RRVM) [14] in a Krylov subspace. 
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2 Rayleigh— Ritz variation method and connected- 
moments polynomial approach 

It is our purpose to solve the eigenvalue equation 

^|^,) = E,|^,), j = 0,l,... (1) 

by means of the RRVM [14] in the Krylov space spanned by | |0j) = H^ \(p)\ . , 

where |0) is a properly chosen vector of the state space. The RRVM makes 
the trial function 

N-l 

\V) = E c, 10,) (2) 

i=o 

orthogonal to the subspace spanned by | |0j) = H^ \(p)\ . ; that is to say: 

{^,\{H~W)\v) = 0,j = 0,l,...,N-l (3) 

where W is one of the RRVM roots and an approximation to the correspond- 
ing eigenvalue. The coefficients Cj are solutions to the secular equations 

N-l 

J2 {H,,~WS,,)c, = 0, j = 0,l,...,iV-l (4) 

2=0 

where the matrix elements Hji = (</)j| H \(f)i) = fii+j+i and Sji = {(f) j\ </)j) = 
/ij+j are given in terms of the moments fij = (0| H^ |0). There are nontrivial 
solutions to the homogeneous system of linear equations (jl]) only for the A^ 
values oi W = Wq, Wi, . . . , Wn-i that are roots of the secular determinant 

\fii+j+i-WiXi+j\f~}^ = 0. (5) 

It is well known that the RRVM approximate eigenvalues approach the 
exact ones from above: W^^^ > W^^+^l > E„, n = 0, 1,...,A^- 1 [14]. 
This statement is not valid if by chance (0| ipi) = because {(f)j\ ipi) = 0, 
J = 0, 1, . . ., and the eigenvector l^/^j) is bypassed by the method because it is 
orthogonal to the Krylov space. Here we assume that the RRVM converges, 
criteria for the convergence of the approach have been discussed elsewhere 
[15]. 



It follows from 
enables us to write 



that (0j| if) = W^ {(f)\ if), which, by virtue of 
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and we conclude that if the coefficients pj satisfy 
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j=0 



0, z = 0, 1, 



N 



then 



N 



0. 



(7) 



(8) 



j=0 

It is clear that we can arbitrarily choose Pq = 1 and that the remaining 
coefficients pj, j = 1,2,. . . ,N are identical to the xat-j+i of Bartashevich 
[13]. 

Since equation ([8]) is satisfied by all the RRVM roots, we conclude that 
it is exactly the characteristic polynomial that results from the secular de- 
terminant dSD- In other words, the CMPA is equivalent to the RRVM which 
explains why the CMPA sequences converge from above [13]. 

In our opinion, solving the well known determinantal equation ([5]) is prob- 
ably easier than solving equations ([7]) and ([8]) of Bartashevich's CMPA. 

If we consider the A^ equations ([7]) plus ([8]) as a whole system of A^ -|- 1 
equations with A^ + 1 unknowns Po,Pi, ■ ■ ■ ,Pn, then we conclude that there 
will be nontrivial solutions provided that 
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fJ'N+1 






(9) 



which is exactly the determinantal equation derived earlier by Bishop et 
al [16]. 

It is interesting to compare the original CMX approaches with the Rayleigh- 
Ritz variation method in the Krylov space discussed above. For concreteness, 
we arbitrarily choose the so-called CMX-LT [3]. Fig[T]shows log \E, 



approx 



E, 



exact 



for the ground state of the dimensionless strongly-anharmonic oscillator H = 
—(P/dx^ + x^. In both cases we chose the trial function {x \(f)) = e~"^ , where 
a = 8, and calculated the "exact" eigenvalue Eq = 1.225820113800492191 by 
means of the Riccati-Pade method that exhibits remarkable rate of conver- 
gence for this kind of models [17]. We appreciate that the RRVM converges 
smoothly from above while the CMX-LT exhibits great oscillation. The same 
situation takes places for other values of a. 

3 Conclusions 

We have shown that the CMPA is equivalent to the RRVM. Consequently 
it applies to all the states and produces sequences of roots that converge 
from above. Numerical difficulties in applications of the CMX have been 
attributed to singularities in the Hankel matrix of connected moments which 
may not be invertible [4] . On the other hand, the Hankel matrix of moments 
that appear in the calculation of the coefficients Xj [13] (or present p/s given 
by ([7])) is always invertible. This is not surprising because one does not 
expect the RRVM to exhibit such kind of singularities. A numerical test on 
a simple anharmonic oscillator suggests that the RRVM may be more reliable 
that the CMX approaches. 
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Figure 1: Rate of convergence measured as log\Eapprox — Eexact\ for the CMX- 
LT (squares) and RRVM (filled circles) in terms of the required number of 
moments M 



